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XX. Coordinate Transformation 

X'n order to facilitate the coincidence calculations, the coor- 
dinates of the ship and the satellite are transformed to the ECO 
system, in which the equatorial plane is the plane of the satellite's 
orbit. The transformation from earth-centered fixed (ECP) to ECO 
coordinates is shown in Figure 1. 

These are the transformation matrices for each step: 

a t -sin 0 

cos. a t 0 
0 0 1 






i&r 
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Step (1) 


Step (2) 
ing node. 


Step (3) 


ECP to ECI. 


tifCi 


Rotate the equatorial plane 
by the Greenwich hour angle 
at time t (0.*.), 



Rotate ECI coordinates by ft, the longitude of the ascend- 



Tilt the equatorial plane by the inclination of the orbit 


Z* 



Figure 1 




cosft 


sinft 0 


(3) X * -sinft cosft 0 


(<l> X 


0 0 
cos i sin i 
-sin i cos i 


The complete transformation is: 


** II I 

(0) X a l A X 


cos (ft - a fc ) 

(-cos i) sin (ft - a^) 
(sin i) sin (ft - a fc ) 


sin (S 1 - a^) 


(cos i) cos (ft - Oj.) sin i 
(-sin i) cos (ft - a fc ) cos i 


The ECO coordinates of the ship are: 


x orb = X cos (^ " + Y sin (ft - a fe ) 


(6) , Y or k = cos i [-X sin (ft - a fc ) + Y cos (ft - a^) ] + Zsin i 
• Z h = sin i [x sin(ft - a.) - Y cos(ft - a.)] + Zcos i 


where X, Y, Z are the ship's ECF coordinates. 
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III. Calculation of Time Intervals 

The ship can be observed when it is in a band about the equator 
in the ECO system. The width of the band is determined by the 

scan pattern of the instrument. The first step in obtaining the times 


the ship enters and leaves the band is to set 2 orfc) 3 + in (6) and 

solve quadratic equations for cc?r and sin a fc (see Appendix 1) . We 
need both sin a fc and cos x { , so that we can put in the proper qua- 
drant. We then use (2) to calculate the U.T. times the ship is at the 


borders of the band in an Interval of approximately 48 hours centered 
about the time of the ship's observation. Depending on the ship's 
latitude, there will be two, three, or four time intervals to be 
searched for coincidences . 



f 
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IV, Existence and Times of Coincidence 


The work of OTIME is nearly complete after the time intervals 
have been calculated. It calls LONGIT, which checks for the exist- 
ence of a coincidence in each interval using the criterion; 

^ on sat ^ fc l^ “ ^ship^l^ (lon sat (t 2 } “ ion ship - 0 ^ 7) 

If the time interval is longer than the period of the satellite, 

LONGIT is restricted to searching the full period closest to the 

time of ship observation. Furthermore, if the interval is longer 

than half a period, it is split into three parts in order to avoid 

a breakdown of (7). 

GETTIM uses a modified regula falsi technique to find the time 
of coincidence. The first step in the procedure uses the normal 
regula falsi technique, The number of iterations to be done is 
passed to the coincidence program from the main program. If conver- 
gence has not been reached after the iterations have been done, the 
second step is to determine in which half of the interval the root 
lies and repeat the first step using that half-interval. The max- 
imum number of times the interval can be halved is also passed from 
the main program. 

If convergence is reached, an empirical correction factor, which 
depends on the month being processed, is applied to the time of coin- 
cidence. The final time of coincidence is 

“ fc c + ‘‘oorrl^o + SorrtV 60 

where t = originally calculated time of coincidence (in hours) 
Sorrl' t corr2 = correction factors. 


( 8 ) 


The correction factors are obtained from plotting <5t, the error in 
the calculated SMMR observation time, against At, the difference 
between the time of the ship observation and the SMMR observation time, 
Points which were observed by SMMR are read from the SMMR TETs for 
three or four days spread over the month being processed, These 
points are fed to OTIME along with several ship observation times. 
These points fit a straight line whose slope, in minutes per hour, is 
^corrl and whose fit-intercept is the negative of t corr2 * 
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V, How to Run the Program 

The following shows how to use OTIME and describes the output. 

Unless otherwise noted, all variables listed below are HEAL* 8. 

A. Input to OTIME 

ALON, ALAT The position of the ship. Units are degrees. 

The range of ALON is (0, 360) , and it is 
measured positive to the west . The range of 
ALAT is (-90, +90). 

TLIM In GETTIM, the interval of coincidence is 

reduced until its length is leas than the 
value of TLIM. A value of ,01 minute is 
recommended. 


NITER (I* 4) 


NHALF (1*4) 


ITSOBS (1*4) 


IDAYYR (1*4) 


LYEAR (1*4) 


Also used in GETTIM, NITER is the number of 
iterations to be done before the time inter- 
val is halved. A value of 3 .is recommended. 

The maximum number of iterations is the pro- 
duct of NHALF and NITER. A value of 17 is 
recommended. 

The hour (U.T.) of the day at which the ship 
made the observation, ITSOBS must be in the 
interval (0,24), 

The day of the year on which the ship made 
the observation. Midnight of January 1, 1978 
is day 0.0. 

The year in which the observation was made, 
LYEAR must have 4 digits, i.e., 1978 not 78. 
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TC0RR1 


Time correction factor a used in GETTXM. 

Units are minutes/hour , TCORR1 is read in 
by the namelist THIS from unit 9. 

TC0RR2 Time correction factor used in GETTXM. Units 

are minutes, TC0RR2 is read by the namelist 
TIME from unit 9, 

When OTXME is called for the first time, it reads in TCORR1 and 
TC0RR2 from unit 9, and from unit 11 it rsads in a set of orbital 
elements and other information which enables it to calculate the 
ECX coordinates of the satellite at any time during the month which 
is being processed. The following cards must be included in the JCL 
for the GO steps 

//GO. FTQ9F0Q1 DD * 
namelist input' 

and 

//GO.FTllFQOl DD DSN-2MMHT . ELTS , aaann , DXSP*SHR 

where aaa » first three letters of the month being processed 
nn * last two digits of the year 

B, Output from OTXME 

TOBS: 4-element time vector containing times of coinci- 

dences measured relative to the time of observation 
(l.e. , « 0). Units are hours (U.T.) . 

IFOBS: 4-element flag vector. If IFOBS(n) » 0, no coin- 

cidence occurred, but if IFOBS (n) =* 1, a coincidence 
occurred, and TOBS (n) contains the time. 

IERF: Error flag (see Sec, C) . 
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C. Error Conditions 


There is on® «rror flag, IERE, which can tak® s®v®rai values s 
Value of TERF Condition 

a. Error conditions sat in OTIHE 

0 Normal 

1 The elements passed to GKEP 

violated one of the follow- 
ing conditions? 

(i) The eccentricity was 
not in the range (0, 1). 

(ii) The semi-major axis 
was not in the range 
(0, 1 x 10 35 ) . 

2 The year in which the 

observation was made does 
not match the year for 
which the orbital elements 
were read in. 


3 


4 

5 

6 
7 


Orbital elements were read 
in for a range of days 
which does not include the 
day the observation was 
made on. 

The hour of the day is out 
of range. 

The latitude of the ship 
is out of range. 

The longitude of the ship 
is out of range. 

The absolute value of a 
computed time is greater 
than 35 hours. 


5 


10 - 


8 The time that the ship 

entered the band is greater 
than the time it left the 
band. 

9 The ship remained in the 

band for more than one 
sidereal day, 

15 One or both of LDAYQ and 

&DAYF, parameters which 
point to the first or last 
days for which orbital 
elements have been calcu- 
lated, are out of range, 
or LDAYO is greater than 
LDAYF. 

16 Orbital elements do not 

exist for the day given by 
IDAYYR 


b. Error conditions set in LONGIT. 

10 Same as 7. 

11 Same as 8. 

12 Same as 9 . 

c. Error conditions set in GETTIM. 

13 The differences in long- 
itude between the ship and 
the satellite have the same 
sign at both endpoints of 
the time interval. 

14 The time of coincidence 

(U.T.) has an absolute 
value greater than 30 
hours . 


NHALF*NITER Convergence was not reached, 

If an error Is encountered, an error message is printed, control 
is returned to the main program, and no ruther processing is done. 


D. Common Block 


The common block CONST appears in all subroutines except GKEP , 
The following table lists the variables contained in CONST. They 
are all computed in OTIME; 


X/ Y, Z 
ALONO 

SREVH 


COSI, SINI 


ECP coordinates of ship. 

Longitude of satellite in ECO system at time 
of ship observation. 

Revolution rate of satellite about the earth 
in radians per hour (U.T.). 

Cosine and sine of orbital inclination. 


Cl, C2 Constants used in calculation of X ^ and 

Y orb ^ see ' and a t ^ see App * 

T Eight element vector containing endpoints of 

time intervals to be searched for coincidences. 


TP I 
BROTH 

GHAO 

TPIE 

SINOM, COSOM 

ERSINI 

PERIOD 

NTIMES 


2 TT . 


Rotation rate of earth, e in radians per 
hour (U.T. ) . 


GHA at time of ship observation. 

2,r/e rot 

Sine and cosine of the longitude of the 
ascending node. 


_e rot ' sin (i) where sin(i) = SINI. 

Period of satellite orbit as computed by 
Kepler's third law. 


The number of time intervals to be searched 
for coincidences is NTIMES divided by 2. 
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E 


Calling Sequence 

CALL OTIME (ALON , ALAT, TOBS, TLIM, IFOBS, IERF, NHALF , NITER, 
ITSOBS, IDAYYR, LYEAR) , 


Appendix 1. Calculation of 

(Al-1) ^ ±z ii ra / s ^ n *•) “ [ x sinQ <- y cosft] cos a fc 

- [x cosfl + y sinQ] sin + z cot i 


(Al-2) 0 *s a cos - B sin a t + C 


where 


A = x sinfi - y cosft 

B = x cosQ + y sinfi 

C = Z coti + Z lim /sin i , 


(Al-3) 


Let sin a t 



cos 



then (Al-2) becomes 


B 


2 


(1 


cos 


2 


a t ) 


2 2 2 
A cos a^_ + 2 AC cosa^ + C 





Appendix 2 . Block Diagram 
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